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ABSTRACT 



For heavy nuclei, very few measurements of the nuclear giant multipole 

resonances have been conducted using inelastic electron scattering. In 

this work, measurements were made to identify the giant multipole resonances 
197 

of Au. Five inelastic scattering spectra were obtained and analyzed for 
beam energies and angles of 75 MeV / 75°, 75 MeV / 90°, 75 MeV / 120°, 

90 MeV / 75°, and 90 MeV / 90°. In each spectrum, giant multipole resonances 
were observed at excitation energies of 7.1, 9.1, 10.9, and 14 MeV. The 
strengths and multipolarities were found to be B(M2) = 15,000e^fm^, 

B(E2) = 8,000e2fm4, B(E2) = 9,000e2fm^, and B(E1) = 90e^fm^ 9 respectively. 

The full widths at half maximum for these resonances were found to be, 
respectively, 3.1, 2.9, 3.1, and 4.9 MeV. 

The resonances found at 10.9 and 14 MeV confirm the only other measure- 
ments made on gold [ 4j . The resonance found at 7.1 MeV agrees with those 
findings concerning the magnetic character of the transition but differs in 
the multipolarity assignment; however, our assignment may also be consid- 
ered tentative since insufficient data were obtained at the backward scat- 
tering angles. The resonance found at 9.1 MeV has not been previously 
reported . 
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I. INTRODUCTION 



The giant dipole (El) resonance has been well known for some time. 

Recent experiments using inelastic electron scattering have revealed the 
existence of other giant resonances in a number of nuclei and in this paper 
we report our investigations about these resonances. 

Many inelastic scattering measurements have been made in the past 15 
years; however, the, giant resonances considered in this paper appear as small 
peaks superimposed on a large continuous spectrum of electrons called the 
elastic radiation tail. Low background and high accelerator beam stabil- 
ity are required to detect these resonances and it is only recently that 
they have been investigated. Studies using Ce, La, Pr [Refs. 1,2] and 
Pb ^Ref. 3j indicate that heavy nuclei tend to have an El resonance at 
about 13 to 15 MeV and an E2 resonance at about 10 to 12 MeV. Other reso- 
nances, some of which were magnetic in character, were observed at excita- 
tion energies ranging from 7 to 9 MeV. Unpublished work at Darmstadt 
[Ref. 4 ] on^^Au has also confirmed this trend: an El resonance was found 
at 10.8 MeV; and one other resonance was located at 7.1 MeV and tentatively 
identified as Ml. Another recent experiment on ^Zr ^Ref. 5 J , found a giant 
resonance (EO or E2) at 28 MeV excitation in addition to the usual E2 and 
El transitions near 15 MeV. 

The investigations reported here were the first attempt to measure 
giant resonances at the NPS LINAC. The procedures and data handling tech- 
niques necessary for the study of giant multipole resonances by inelastic 
electron scatterings have been developed and refined. 

Experimental values for the inelastic scattering form factors were 
obtained using incident energies of 75 and 90 MeV and scattering angles of 
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75, 90, 120, and 150 degrees. The results of this work agree with other 
reported works ^Ref. in that the previously reported resonances were 
observed. A previously unreported resonance was found at 9.1 MeV with 
B(E2) = 8,000e^fm^. Additionally, preliminary indications for the exis- 

i 

tence of an E2 resonance in the vicinity of 23 MeV were obtained. 

This thesis is organized into five sections and one appendix. In 
Section II a brief review of inelastic scattering and giant resonance 
theory is presented. Section III briefly describes the NPS LINAC and the 
experimental procedure used by the authors. Section IV discusses the meth- 
ods employed for data reduction. Topics include radiation corrections, 
calculations of the elastic radiation tails, procedures used in the loca- 
tion and characterization of the giant resonances, and error analysis. 
Experimental results and conclusions are presented in Section V. The cal- 
culations involved in this work are quite complex and cannot reasonably be 
performed without the use of various sohpisticated computer programs. The 
appendix contains the program listings and detailed explanations concerning 
their use. 
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II. THEORY 



A. INELASTIC ELECTRON SCATTERING 

When a relativistic, mono-energetic electron beam is incident on a 
thin target most electrons pass through, but some collide elastically with 
target nuclei and a few collide inelastically , thereby exciting nuclei to 
higher allowed energy states. 

The energy distribution of the scattered electrons is determined ex- 
perimentally at various scattering angles. Most of the electrons observed 
at a given energy, E^ , less than the incident energy, E^, are due to elas- 
tic collisions with nuclei accompanied by energy loss in passing through 
the target. Energy loss mechanisms include: (1) ionization losses, 

(2) emission of a photon during the scattering event, and (3) bremsstrah- 
lung while passing through the target both before and after the elastic 
scattering event. This "elastic tail," and any other "background” affect- 
ing the observed spectrum, must be calculated and subtracted from the ex- 
perimental energy spectrum to obtain the inelastic spectrum. The resulting 
unfolded inelastic spectrum consists of the inelastic peaks and their 
associated tails. Each inelastic peak corresponds to an excited state of 
the nucleus. 

Inelastic collisions are classified into three types, according to the 
mechanism of interaction with the nucleus: 

1 . Surface Nucleon Excitation 

One or more surface nucleons are excited with excitation energies 
ranging from several KeV to a few MeV. The higher energies tend to be 
associated with collective excitations and have larger characteristic 
resonance widths. 
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2. Giant Resonances 



Many nucleons participate in oscillations causing the giant reso- 
nances. These resonances have very large natural widths and are found 

r 

typically in the region of ten to 20 MeV excitation for heavy nuclei. This 
large resonance width occurs because the excitation energy is large enough 
that the excited nucleus can decay by particle emission. The short life- 
times associated with particle emission result in the large widths because 
of the Heisenberg uncertainty relationship between energy and time. 

3 . Quasi-Elastic Resonances 

In this process a nucleon is ejected from the nucleus during the 
collision. The incident electron interacts with essentially only one 
nucleon. These "quasi-elastic" collisions produce a very wide peak in the 
spectrum at an excitation energy of a few hundred MeV for heavy nuclei. 

B. GIANT RESONANCES 

According to the Generalized Goldhaber-Teller Model ^Ref. 6j , nuclei 
are capable of "collective oscillations" in which all the nucleons oscil- 
late inside the nucleus. There are four separate fluids in motion: neutrons 

with spin up, protons with spin up, neutrons with spin down, and protons 

with spin down. Each fluid is characterized by the spin and the isospin 

of its particles. 

The motion of the four fluids, two of which are 180° out of phase with 

the other two, can be characterized by two numbers, S and T. If protons 

and neutrons move in phase, T=0; if they move 180° out of phase, T=l. If 
particles with spin up are in phase with those with spin down, S=0; if they 
are 180° out of phase, S=l. A transition to an S=1 state changes the mag- 
netization and is therefore called a magnetic transition. A transition to 
an S=0 state is called an electric transition. 
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The electric and magnetic fields of nuclear excited states are described 
in terms of multipole expansions. A transition is labeled according to the 
term in the multipole expansion which is predominantly affected. For ex- 
ample, a transition from the ground state which increases the magnetic (or 
electric) dipole term in the electromagnetic multipole expansion is labeled 
Ml (or El). If the electric (or magnetic) quadrupole term is involved, the 
transition is labeled E2 (or M2). An electric monopole transition is labeled 
EO. There are no magnetic monopole transitions. In general, L units of 
angular momentum are required for an EL or ML transition from the ground 
state. 

C. NUCLEAR FORM FACTORS 

1 . Development of the Form Factor 

In the plane wave solution to the non-relativistic wave equation 
we find that the differential elastic scattering cross section can be ex- 
pressed as the magnitude of the square of some function of scattering angle, 
f(0), where 

W-^/lXOe'^d’r, 

_2m ^ 

and U(r) = ^2 V(r), where V(r) is the potential of the electron in the 
field of the nucleus, and K is the wave vector of the incident wave. 
Replacing U(r) by the related electrostatic potential, <p( r) , through the 
relation V(r) = q 0(r) , we have 

f («)=♦$ ff e‘* -1 V>dV. 

We can relate 0(r) to the charge distribution p( r) through Poisson’s equa- 
tion; taking the fourier transforms and equating integrands leads to the 
result 
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f(«)= -|f#- 



with 

i>(K)= {Tiffpine*" dV = (gV^ZeF^). 

Hence , 

f(ffl = F(S). 

We can now solve for the differential cross section to obtain 

do__ 4m 2 eV|p ( K)| 2 
dn (h k) 

For a point charge it can be shown that 5 S ° we ^ ave 




We call F(k) the form factor for the charge distribution with 

F(«) = Yefi p(?)e’‘ K *^d 3 r . m-D 



For inelastic scattering, the charge distribution changes during the 
transition. In this case, p(r) is replaced by the four-vector Dirac cur- 
rent operator, , associated with the transition involved. 

2. Mott Cross Section 

The differential cross section per unit solid angle for elastic 

scattering of a relativistic electron by a point nucleus of charge Z was 

shown by Mott [Ref. 7] to be: 

Ida ] _ Z 2 e 4 COS*(9/2 p = /da\ p 

\dn/ pt 4E- sin 4 0/2 * ldn] m ott* 7 

where R is the recoil factor defined by the relation 
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1 



R* 1 + 2Ef 

+ Mc 2 sin 2 0/2 



Using the concept of the form factor we are able to generalize this result 
to finite charge distributions to obtain 



M»\ _ M®.] . I 

Idnj [dill rnoti 

For heavy nuclei where E^< 100 MeV, R = 1 ; 




therefore 




mott 




(H-2) 



If this expression is used with the form factor given by equation II-l, 
we have the Born approximation for elastic scattering. For large Z and low 
energy this approximation is inaccurate and the cross section must be cal- 
culated by partial wave analysis [Ref. 8]. 

3. Electron Scattering from a Proton 

The electromagnetic interaction Hamiltonian between the electron 
and a point proton is expressed in four-vector notation as 
Calculating in terms of the Dirac theory operators [Kef. 9] gives 

V- 



Solving the Dirac wave equation yields the point Dirac proton solution: 

ck = /ck] . R /, + jl^tanVa) =/ds.] •R-|F,(q 5 )f 

dn Idn/mott r idft/mott '' 1 ' 



The first term in the parentheses corresponds to the electric interaction 

and the remaining term corresponds to the magnetic interaction. General- 

2 

izing to finite charge distributions and adding another term, F£(q )> to 
account for the anomalous magnetic moment distribution we obtain 
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da _/da \ 

dQ \dn/ mott 



R(|«i 



* + q 



4 m 2 c 2 



+ 



q 2 



4 m 2 c 2 



tan^lF^/iF J ). 



The first two terms account for the electric interactions and the last term 
accounts for the magnetic scattering; q 7 is the square of the four-momentum 
transfer and /z is the proton magnetic moment in nuclear magnetons. 

This result, although it does not apply directly to these investigations, 
does show the presence of two form factors, due to scattering from the charge 
and magnetic moment of the proton. 

4 . Plane Wave Born Approximation (PWBA) 

The Born approximation is used to calculate inelastic scattering 
cross sections based on an assumed nuclear transition charge distribution. 

The PWBA uses first order perturbation theory to calculate cross sections 
for interactions between an electron and the electromagnetic field of the 
nucleus. A multipole expansion is used to describe the fields. The tar- 
get nuclei are assumed to recoil non-relativistically . The PWBA assumes 
plane wave solutions for the incident and scattered electron. 

For E^, E£»Ex, Isabelle [Ref. lo] has shown that the differential 
inelastic cross section per unit solid angle is given by 

M; ft) mott fc q ' L) + ( * +,an2 ‘ /2)F ’ (q ' L >] , 



and 



'do\ = tda\ 

,dn/ ML \dn/ mott 



[( 2 + tan 2 0/2) F^cq!*] ^ 



where F^, and are the form factors for longitudinal and transverse photon 
absorption. For elastic scattering, the PWBA results are equivalent to 
equations II- 1 and II-2. 
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5 . Distorted Wave Born Approximation (DWBA) 



For heavy nuclei, the coulomb field surrounding the nucleus dis- 
torts the electron wave function sufficiently that a plane wave solution 
is inadequate. The DWBA assumes phase-shifted spherical waves for the in- 
cident and scattered electron wave functions. Unlike the plane wave Born 
approximation, a form for the charge distribution function must be assumed 
before calculations are begun. 

For elastic scattering, the cross section may be calculated using a 
Fermi charge distribution. However, this is not appropriate for inelastic 
scattering, in which the outer nucleons only would be expected to take part 
in a nuclear transition. A good approximation was suggested by Tassie 
[Ref. ll] , in which the transition charge distribution is assumed to be 
the derivative of the Fermi function. Computer programs for calculation 
of the differential cross sections (^) ei and (^) ML have been devised by 
Tuan, Wright, and Onley []Ref. 12]]. 

D. IDENTIFICATION OF EXCITATION MECHANISMS 



The form factor squared for an inelastic transition is defined by 



H -M 

ldn/m.„u Idn/mott 






Using the DWBA, the theoretical variation of F^ with scattering angle can 
be calculated and plotted for each EL and ML at any incident energy. 

In actuality, one needs to consider only L~0 through about L=3 since 
higher terms are expected to contribute a neglible amount to the inelastic 
scattering. Each observed peak of the inelastic spectrum can then be 
identified by plotting the experimental form factor variation with angle 
at each incident energy used in the experiment, and comparing with the 
theoretical curves. 
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III. EXPERIMENTAL EQUIPMENT AND PROCEDURE 



A. NAVAL POSTGRADUATE SCHOOL LINEAR ACCELERATOR 

The NPS 120 MeV LINAC was described in detail by Barnett and Cunneen 
[Ref. 13] . It is an electron accelerator similar in design to the Stanford 
University Mark III linear accelerator described by Chodorow et . al . 

^Ref. 14j . The 30 foot NPS LINAC is powered by three klystrons. Each 
klystron can yield up to 22 megawatts of peak power. The machine is pulsed 
at a rate of 60 pulses per second and each beam pulse is approximately one 
microsecond in duration. The maximum average beam current is approximately 
one microampere. 

Beam deflection magnets and energy defining slits are used to select 
the desired beam energy, which may range from 30 to 105 MeV. The field of 
the beam deflection magnet is measured by a nuclear magnetic resonance 
(NMR) probe. Quadrupole magnets focus the beam to a diameter of approxi- 
mately two millimeters at the target. 

B. SPECTROMETER AND COUNTING SYSTEM 

The NPS LINAC Spectrometer was described in detail by Oberdier jjlef. 15] . 
It is a movable, 16-inch radius, double-focusing magnetic spectrometer 
capable of measuring electrons scattered at angles from 45° to 150° in 
15 degree stops. 

The field of the spectrometer magnet is varied to measure the energy 
distribution of the scattered electrons. Electrons with energies within 
— 1 . 5 per cent of the spectrometer energy setting are detected by one of 
ten NE-102 plastic scintillation counters located on the focal plane of 
the spectrometer magnet. Each detector serves as the source of a counting 
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channel. Each counter has a width which allows detection in an energy spread 
of 0.3 per cent of the central spectrometer energy. Two large detectors, 
each of which spans the same field as the ten smaller detectors, serve as 
dual coincidence backing counters. A charged particle passing through a 
front detector and both backing detectors will register as one count on a 
Jorway Model 1836 Dual Channel Scaler located in the control room. The 
counter gating circuit is synchronized with the accelerator gun pulse so 
that those triple coincidences that occur between pulses are not counted. 

The scalers are capable of recording periodic pulses up to 150 MHz. 

However, since the beam is pulsed, the count rate was limited to 50 counts 

per second to avoid significant count' rate losses. A Faraday cup is used 

to calibrate the charge collected on the Seconday Emission Monitor (SEM) . 

197 

Three pure Au foil targets were used with thicknesses of two mils, 
five mils, and ten mils. The targets were placed in an evacuated chamber 
with 1.0 rail Havar windows positioned to allow the incident electrons to 
enter and the scattered electrons to exit. Since the scattering cross sec- 
tions decrease at the backward angles, five and ten-mil targets were used 
in order to increase the count rates and thereby reduce the effect of sta- 
tistical fluctuations. The two-mil target is recommended for future in- 
vestigations to maintain adequate resolution (the ten-mil 

target gave an elastic peak width of 19.5 per mil while the two-mil target 
gave an elastic peak width of 6.7 per mil.) 

C. DATA ACCUMULATION 

During each experimental "run" the beam energy and scattering angle 
were held constant. The number of counts in each of the ten counting chan- 
nels and the integration capacitance and voltage were recorded automatically 
at each spectrometer energy. The spectrometer energy was decreased in steps of 
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0.1 MeV. The integration capacitance and voltage were chosen to maintain a 
high count rate not to exceed 50 counts per second. The beam current was 
decreased when measuring the elastic peak to avoid an excessive count rate. 
The spectrometer energy was reduced from 100 MeV to about two MeV above the 
elastic peak. Elastic scattering data were taken at 0.1 MeV steps to about 
four MeV below the elastic peak. When the elastic data were accumulated, 
the beam was retuned to give maximum current allowed by count rate restric- 
tions, and the spectrometer energy was set back several MeV to overlap the 
elastic data. In order to insure that the magnet always followed the same 
hysteresis loop, the spectrometer magnet was returned to a field value cor- 
responding to 100 MeV and slowly decreased to the desired energy. It was 
not necessary to consider relative counter efficiencies in this work. The 
KIK code [Ref. 16] automatically corrects the energy of each counter to 
correspond to the true energy being detected. 
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IV. DATA ANALYSIS 



A. DATA REDUCTION 

The data were collected on a Model 33 teletype on paper and punched 
tape simultaneously. The punched tape was transcribed onto computer data 
cards compatible for use on the IBM 360-C. The organization of the data 
cards is described in Appendix A-2. Fortran program DATA LIST was used to 
print out the transcribed data cards for verification against the teletype 
"hard copy." After verification and correction, if required, the data 
cards were processed by fortran program KIK. This computer program assimi- 
lates the ten spectra (one spectrum from each of the ten spectrometer chan- 
nels) into one total spectrum, computes the statistical error, and outputs 
the results in both hard copy and punched card form. The punched cards 
from KIK are used directly by fortran program SPECTRUM to draw the total 
spectrum as a point plot on the CALCOMP plotters at the NPS Computer Center. 
These spectra are displayed as Figures IV-i through IV-5 . 

* The spectrum cards from KIK were then processed by fortran program NAW 
[kef. 16] . This program best fits an analytic expression to the elastic 
peak data and computes the area under it - the integration is terminated at 
AE = 2r s ,whereQ is the full-width of the elastic peak measured at half- 
maximum. NAW then makes the appropriate radiation corrections to the area, 
as detailed in Section IV-B. Experience has indicated that the best value 
for parameter L^ is L^=3 (see Appendix A-4). 

Program NAW then determines an analytic expression to represent the 
background radiation due to the elastic tail in order to make the requisite 
subtractions. The necessity for these subtractions will be discussed later 
in Section IV-C . Experience of experimenters at Institut Fur Technische 
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Kernphysik at Darmstadt led to the choice of parameter 1^=3 and L^=5 
resulting in the "best" representation of the radiation tail. The radi- 
ation tail and any constant background (measured above the elastic peak) 
are subtracted from the data, and program NAW then determines a best fit 
analytic expression for the inelastic resonance peaks. The technique used 
was to fix the positions of the resonance peaks based on viewing the various 
spectra, select peak widths (full-width at half -maximum) based on other 
published experimental results, and then dynamically adjust these parameters 
until one best fit for all experimental spectra was obtained. Program NAW 
displays these results both graphically and in tabular form. Furthermore, 
it integrates the analytic expression for the inelastic resonances to deter- 
mine the areas under the peaks and displays the ratio of the area of the 
inelastic peak to the area of the elastic peak for each resonance. 

Using the incident electron energies calculated by NAW, fortran program 
MARPHASE uses the Rawitscher-Fischer Partial Wave Code [Ref. 8] to calcu- 
late the elastic scattering parameters; the results are displayed in 
Table IV- 1. Since the differential scattering cross sections are propor- 
tional to the area under the resonance peaks one can easily determine the 
inelastic scattering cross sections, or more usefully, the square of the 

o 

inelastic scattering form factor, F : 

p 2 = (<f<r/dJi)inet Ac (dr/d.st)efa 5 

(d <r/ d -ft) roott ^ ^ (d<T /cl • 

The inelastic form factors squared may now be graphed as a function of the 
scattering angle and parametized by incident energy. To identify the ex- 
citation involved, one may compare these angular plots with the theoretical 
values as determined by a DWBA code. Fortran program DUKE [kef. 12] uses 
a DWBA code to calculate F for both electric and magnetic transitions. 
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TABLE IV- 1 



ELASTIC SCATTERING PARAMETERS FOR 197 AU 



0 


E t - (MeV) 


da/dfl 


(cWdn) Mi , t 


F 2 


75° 


75.15 


2.52x10"* 


2.62 


9.59xl0~ 2 


75° 


90.19 


7. 83x1 0" 2 


1.82 


4.30x10-2 


90° 


75.20 


5.62xl0- 2 


1>14 -1 
7.95x10 1 


4.91x10-2 


90° 


90.19 


2. 39xl0“ 2 


3.00x10-2 


120° 


75.06 


8.75xl0- 3 


2.55xl0 _1 


3.43x10-2 



Differential cross sections are in units of fib/sr. 



A transition strength, B, of unity is assumed for all transitions in these 

calculations. While the output labeled as F0RM is the theoretical value 

9 2 
for for electric transitions, the corresponding output is F multiplied 

by the factor V^/V^ (labeled QVLT on the output) for magnetic transitions. 

Once the theoretical values for F have been obtained the experimental 

values for the transition strength can be calculated by simply averaging 

the values, i.e. 

Fortran program BCALC computes B and then multiplies the theoretical values 

2 - 2 

of F by B. These are graphed, and the experimental values of F are then 

plotted on the same graph. This procedure is done with each identified 
resonance for all reasonably possible transitions. The transition type 
identification is made by determining which graph exhibits the best fit of 
the experimental points to the theoretical curves and also yields a reason- 
able value for B. To avoid subjective, and thereby possibly prejudiced, 
evaluations of the best fit, program BCALC computes the normalized quanti- 
ties', VAR, which are compared to determine the best fit - the smaller the 

value of VAR, the better the fit. x 
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VAR =2('S^y^). 



€tp VQ^'i 

These values may be found in Table V-2. 



B. RADIATION CORRECTION CALCULATIONS 

There are processes occurring during elastic scattering that result in 
low energy electrons which constitute the so-called elastic radiation tail. 
Since it is not practicable to integrate the electron spectrum under the 
elastic peak to zero energy, it is conventional to apply three correction 
factors. It is sufficient to correct for such processes as follows: 



feb ffiw. . 



where Kg, Kg, and K^ are, respectively, the Schwinger, Bremsstrahlung, and 
Electron-Electron Scattering corrections. 

1 . Schwinger Correction 

The Schwinger term is by far the most significant of the three 
corrections being discussed. This term corrects for virtual photon emission 
and reabsorption and real photon emission occurring in the field of the tar- 
get nucleus - both of which result in large-angle scattering of the incom- 
ing electron, i.e., during the scattering event. Schwinger's original 
work [Refs. 17, 18] gives the relation 



K s - 






where a is the Sommerfield fine structure constant, E^ is the most probable 
energy of the outcoming electron, AE is the width of the integration (in 
program NAW AE is set at twice the elastic peak full-width at half-maximum), 
and f(0) is the Spence Function. Many investigations have since improved 
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upon this formulation, notably the work of Tsai [^Ref. 19] and that of 
Meister and Yennie [jlef. 20] . 

2 . Bremsstrahlung Correction 

This term corrects for the energy radiated via bremsstrahlung during 
small-angle deflections caused by the electron’s interactions with the fields 
of nuclei other than the elastic scattering center both before and after the 
elastic scattering event* This energy loss mechanism was first studied by 
Bethe and Heitler [Ref. 2lj , and the correction factor for this mechanism 
was given by Hofstadter [Ref. 22] as 



t< 8 = exp 



t 

XJn 2 







where X 0 is the radiation length and t is the target thickness, both in 
g/cm^. Program NAW uses the formalism of Barber, trt. al . ^Ref. 23^j to give 



= 



= expj— -t-f—L— % 
P |2cos 0/2^X o (E 1 -E f ) E i 



(E-2 + Ep 2) 2 1 

E j .E f 3 91n(183exp{-lnz Vi }) 



3. Electron-Electron Scattering 

This term, sometimes called the Moller term, is almost negligible. 
The term corrects for electron-electron scattering in the target that occurs 
before or after the elastic scattering event. The effect is to produce low- 
er energy electrons, almost in the original direction, similar to the effect 
of the small-angle bremsstrahlung. By Barber’s formalism we have 

~1 



Km = 



exp 




2E_l + 1 



_ , „ a. 1 , 1 

Ef E f (E i ~E f ) (E^+l ) ^ + TET^F + TeT+TF 



C. RADIATION TAIL CALCULATIONS 

When calculating the radiation corrections we were concerned only with 
determining how many elastically scattered electrons were degraded to an 
energy beyond the integration cut-off, Ae. The calculation of the radiation 
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tail is concerned with the energy distribution of electrons lost from the 
elastic peak. Program NAW has the option of using either Schiff's method 



The latter was used in this experiment resulting in an approximation to the 
true radiation tail - while the approximation is crude, it does have the 
proper slope. NAW then best fits a f our th~degree polynomial to the cal- 
culated values (output in program NAW as SSW) . The resulting expression 
for SSW is then fit to the raw data spectrum after any constant background 
counts have been subtracted. Of the 11 fit options available, the one used 
in this work was 



where: UNT is the portion of the measured counts//iC attributed to the 

elastic radiation tail; A, B, and C are numerical fitting parameters deter- 
mined by NAW; and X = E-Eg (Eg is a constant determined by NAW). The radi- 
ation tail obtained is then subtracted from the inelastic data to yield 
the unfolded spectrum. 

D. ERROR ANALYSIS 

1 . Systematic Errors 

In this experiment, the ratio of the inelastic to the elastic cross 
section is computed. Since the inelastic and the elastic data were always 
taken during the same experimental run, most systematic errors cancelled. 
These would include errors due to variations in solid angle, scattering 
angle, SEM efficiency, target thickness, target purity, and absolute spec- 
trometer efficiency. On several occasions after the inelastic data points 
were collected, the elastic peak points were repeated to verify that the 
number of counts per microcoulomb remained constant, within statistical 
variations, during a run. 




24 or the more recent method 




UNT = A* SSW + B + C-X 
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2. Random Errors 



The statistical fluctuations of the number of electrons scattered 
for a given number of incident electrons obeys Poisson statistics, in which 
the standard deviation is given by S^CE) = ^ N (E) , where N(E) is the total 
number of counts at scattered energy E. The fractional error is therefore 
given by: 

J N(E)' 

Fractional error = 7 —-. 

N(E) 

Program KIK computes the energy, E, for each counter at each spectrometer 
setting, sums the ten spectra to yield one total spectrum, and tabulates 
the total counts per microcoulomb and the associated fractional error. 

This error varied from one percent to three percent. The effect of this 
error on the final form factors is calculated by NAW as follows. Let the 
total counts detected in an energy region encompassing the width of an in- 
elastic peak be and the counts attributed to the inelastic peak be 
(thus - Nj is attibuted to radiation tail and background in this region) . 
Then the fractional error in the cross section is given by \J^t' / n I* ^ ue 
to the varying strength of each transition and the fact that the total num- 
ber of counts was lower for certain runs, these errors varied greatly in 
size. They are tabulated in Table V-l. 

Other random errors include changes in the beam current, variation in 
SEM efficiency during a run, energy shifts of the incident beam, leakage 
current in the beam current integrators, error in measurement of the back- 
ground, uncertainties in the position of the peaks, and uncertainties in 
the widths of the peaks. 

Variations in the beam current affect the values obtained for the 
cross sections since count rate losses increase with current. Although 
these losses were initially believed small for count rates up to 50 counts 
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/ 



per second, current studies indicate that this is probably not so. The 
magnitude of this error is, as yet, unknown, but may well contribute as 
much as an additional 25 percent error to the experimental values of the 
form factors. 

The possibility of variations of the SEM efficiency during a run was 
not investigated in this experiment. Measurements by Berard and Traverso 



Variations in line voltage cause fluctuations, not only in beam energy, 
but in spectrometer and beam deflection magnet fields as well. The effect 
is to shift some electrons from one counter to another. However, the elec- 
tronic regulation of the spectrometer and deflection magnets reduce this 
error to a negligible level. Leakage currents in the beam integration cir- 
cuit were measured by Topping J^Ref. 27^ and were found to be negligible 



The background contributes less than 1 percent of the area of the elas- 
tic peak, so that the fractional error due to the background subtraction 
is negligible - particularly when compared to the uncertainty involved in 
the subtraction of the radiation tail. 



to assess. The experimental form factors calculated by program NAW depend 
on the choice of location and width of each of the inelastic peaks. It is 
estimated that these uncertainties contribute a 20 percent error to the 
experimental values of the form factors. 

The total fractional error for any given calculation of an experimental 
form factor can now be obtained as follows: 




show that the SEM efficiency remains constant within 0.09 percent. 



also. 



The error involved in fitting inelastic peaks to the data is difficult 




a\\ ttvrct | 
trrer 
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where is the fractional error due to the i’th mechanism. Since three 

factors dominate the others, the fractional error can be approximated by 
f o 2 9" 1 

S = J(.25) + (.2) + , where is the fractional error due to Poisson 

counting statistics (found in Table V-l). In most cases, the total esti- 
mated error was approximately 40 percent. 
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V. EXPERIMENTAL RESULTS AND CONCLUSIONS 



A. EXPERIMENTAL RESULTS 

The results of five experimental runs are reported here. Three other 
runs were conducted but are not reported. Two runs conducted at a scatter- 
ing angle of 150° had counts that were so small that any resonances present 
are masked by the statistical fluctuations. One run conducted at an inci- 
dent energy of 90 MeV and a scattering angle of 120° using a ten-mil target 
had insufficient resolutions to extract any structure from the resultant 
spectrum. Prior to presenting the experimental results, a brief description 
of the procedure by which program NAW was used follows. 

The initial positions of the resonances were estimated by viewing the 
experimental spectra, resulting in the values of 7.4, 8.8, 10.6, 13.8 and 
14.2, and 23.0 MeV. Since the fit by NAW for the radiation tail is only 
considered valid for a range of about 20 MeV and the theoretical values for 
the radiation tail are known to become progressively inaccurate with in- 
creasing excitation energy, the 23.0 MeV resonance was not investigated in 
detail (preliminary results indicate that this resonance is probably an 
E2 transition). The remaining four resonances were given to program NAW, 
with widths of 3.0, 2.5, 3.0, and 4.75 MeV respectively, and the program 
was allowed to vary the positions of the last two resonances. Each run 
resulted in somewhat different values for the positions of the resonances, 
and average positions were determined. The averaged values were then en- 
tered into program NAW and, keeping all four resonances fixed in energy, 
the widths were all allowed to be determined by NAW. At this stage, for 
each resonance, the average value of the width for all the runs was calcu- 
lated, Entering these averaged values for the widths into NAW, the 7.4 
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and 8.8 MeV positions were allowed to vary and again the averaging process 
was followed. The new positions were given to NAW, the widths were once 
more allowed to vary, and the results were again averaged. The final re- 
sults of this procedure were resonances at 7.14, 9.12, 10.9, and 13.8/14.2 
MeV with respective widths of 3.13, 2.89, 3.06, and 4.89 MeV. The peak in 
the vicinity of 14 MeV was found to be incident energy dependent. At an 
incident energy of 75 MeV the peak maximum was located at 13.8 MeV, and at 
90 MeV the peak maximum was located at 14.2 MeV. Hereafter this resonance 
will be referred to as the 14 MeV peak. 

Program NAW determined an analytic expression to best fit the data, 
resulting in a smoothed spectrum. In also determines an analytic expres- 
sion for the radiation tail; these functions are displayed graphically in 
Figures V-l through V-5 . The radiation tail and any constant background 
were then subtracted from the smoothed spectrum, and analytic expressions 
for the four inelastic resonances were determined assuming a Breit-Wigner 
form. The areas under the inelastic resonances were then computed and the 
ratio of each of these areas to the area under the elastic peak was calcu- 
lated along with the fractional error due to counting statistics (Poisson) . 
Multiplying A^/A e by the square of the elastic scattering form factor (cal- 
culated by program MARPHASE and displayed in Table IV-1) gives the square 
of the inelastic form factor. The results of these calculations are dis- 
played in Table V-l. 

2 

Using the experimental values of F at the different angles and ener- 
gies and the theoretical values as calculated by DUKE as input data, program 
BCALC determines the transition strengths, B, and the values of VAR for 
the various transitions. It also prints out (for subsequent plotting) 

o 

experimental and theoretical values of F based on the calculated value 
of B. The numerical results are displayed in Table V-2 and the graphs are 

reproduced as Figures V-6 through V-9. 
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TABLE V-l 



EXPERIMENTAL INELASTIC SCATTERING PARAMETERS FOR 197 AU 



E x (MeV) 


t (mil) 


0 


E Q (MeV) 


Aj/A e 


Percent 
Er ror 


F 2 




2 


75° 


75.04 


2 . 24xl0 -2 


24 


2.14xl0“ 3 




2 


90° 


75.07 


2 . 41xl0 -2 


30 


1 . 18xlO~ 3 


7.14 


5 


120° 


74.82 


3.45xl0~ 2 


12 


1 . 18x1 0 -3 




2 


75 ° 


90.07 


3.54x10 2 


17 


1 . 52xl0 -3 




2 


90° 


90.05 


2.03xl0 -2 


24 


6. lOxlO -4 


/ 


2 


75° 


75.04 


7 . 07xl0 -3 


18 


6. 78xl0 -4 




2 


90° 


75.07 


7 . 05xl0 -3 


24 


3.46xlO“ 4 


9.12 


5 


120° 


74.82 


4.40xl0" 3 


41 


1 . 51xl0 -4 




2 


75° 


90.07 


1.20xl0- 2 


11 


5 . 16xl0 -4 




2 


90° 


90.05 


6. 93xl0 -3 


22 


2 . 08xl0 -4 




2 


75° 


75.04 


1.06x10-2 


9 


1 . 02xl0 -3 




2 


90° 


75.07 


l.lOxlO -2 


13 


5 . 38xl0 -4 


10.90 


5 


120° 


74.82 


2 . 08xl0 -3 


64 


7 . 14xl0 -5 




2 


75° 


90.07 


1.37x10-2 


9 


5 . 90xl0 -4 




2 


90° 


90.05 


8. 79xl0 -3 


16 


2. 64x1 O -4 




2 


75° 


75.04 


1.56x10-2 


7 


1 . 50xl0 -3 




2 


90° 


75.07 


1. 32x10-2 


14 


6 . 49xlO -4 


14.00 


5 


120° 


74.83 


5. 69x1 0 -4 


298 


1.95xl0 -5 




2 


75° 


90.07 


1.59x10-2 


10 


6 . 84xl0 -4 




2 


90° 


90.05 


l.llxlO -2 


15 


3. 32xlO -4 



TABLE V-2 

MULTIPOLARITY PARAMETERS 



m 

X 

/*> 

3: 

a> 

< 

'w' 


Type 


B(e‘fm lL ) 


VAR 


7.14 


Ml 


163.9 


0.956 


7.14 


M2 


14,895.9 


0.511 


7.14 


E2 


31,655.0 


0.880 


9.12 


Ml 


48.4 


1.574 


9.12 


E2 


7,819.7 


0.125 


10.90 


El 


91.2 


1.105 


iO . 90 


E2 


8,822.9 


0.582 


14.00 


El 


89.9 


1.502 


14.00 


E2 


10,527.6 


1.390 
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B. CONCLUSIONS 



Transition assignments were made as discussed in Section IV-A. The 
7.14 MeV resonance appears to be an M2 transition, but this assignment can- 
not be made conclusively due to lack of scattering data at the background 
angles. Future investigations should concentrate on obtaining data at the 
backscattering angles, since the cross sections for magnetic transitions 
increase at these angles. It may also be useful to obtain data at the small 
forward angles, in order to move away from the peaks in the form factors. 

The heretofore unreported resonance at 9.12 MeV is concluded to be an E2 
transition. The assignment of the 10.9 MeV resonance as an E2 is in agree- 
ment with previous unpublished work at Darmstadt [kef. 4]. The results of 
our work indicate that the resonance collectively referred to as 14 MeV 
shifts with incident energy. This shift has been explained by Pitthan 
[Ref. 2^ , who assumed that the broad resonance in the cross section is cal- 
culated from a broad nuclear excitation strength function. This strength 
function varies primarily with excitation energy, but also varies with in- 
cident beam energy and scattering angle. Thus the maximum of the inelastic 
peak may shift when beam energy is changed. The results of this work, 



incident energies) indicate this resonance is an El giant resonance. 

The results described above are recapitulated in Table V-3. All pre- 
viously reported j^Ref.4j resonances appeared in our spectra. The observed 
strengths were not in exact agreement with the previous results, but dis- 
agreement can be explained by observed deficiencies in the counting system. 
Furthermore, two additional resonances were observed: 9.1 MeV, with 

B(E2) = 8,000e^fm^; and 23 MeV which is probably an E2 transition. The 
limitations of the determination of the elastic radition tail, and insuf- 
ficient experimental data precluded detailed investigations of the latter 



coupled with the results of the Darmstadt group (conducted at lower 
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resonance. We have demonstrated that giant resonances can be measured at 
the NPS LINAC Facility. The various experimental techniques developed, the 
problem areas identified and rectified where possible, and the computer 
data handling procedures developed will enable future investigations of 
giant resonances by the LINAC group to proceed much more rapidly and accu- 
rately. Future investigators should spend a fair amount of time developing 

new techniques for fitting the elastic radition tail, because its effect 

2 

on the resulting values for F is paramount. 



TABLE V-3 

IDENTIFIED GIANT MULTIPOLE RESONANCES IN 197 AU 



E x (MeV) 


Q (MeV) 


Transit ion 
Type 




7.1 


3.1 


M2 


15,000 


9.1 


2.9 


E2 


8,000 


10.9 


3.1 


E2 


9,000 


14.0 


4.9 


El 


90 
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Figure V - 1 
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Figure V-2 
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Figure V-3 
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Figure V-5 
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Figure V-8 
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APPENDIX A 



COMPUTER PROGRAMS 



Seven FORTRAN programs were used to perform the necessary calculations 
for this work. The computer programs are reproduced on the following pages. 
Each program listing will be prefaced by a brief explanation of what it does, 
and a parameter list with explanatory remarks (program options not directly 
applicable to this and related works will not be listed). Additionally, if 
any special computer control cards are required to execute the program at 
the NPS Computer Facility, a sample FORTRAN deck set-up will be included. 

A. DATA LIST 

DATA LIST prints the raw data cards in the teletypewriter format used 
by the NPS LINAC Data Acquisition System, segregating by integral data sets, 
so as to aid in verifying the accuracy maintained during the processing from 
punched-tape output to punched-card output. This program requires only one 
parameter card: 

TEST, TEXT (A4,19A4) 

TEST The first four characters on the first data card. 

Usually a blank followed by three digits. 

TEXT Any alphanumeric message to be printed on the first 

page of the data listing. 

B. KIK 

KIK can process up to 499 data sets, each set consisting of the scaler 
readings for each of the ten spectrometer channels of the NPS LINAC Count- 
ing System. KIK compiles the ten separate spectra into one total spectrum 
and calculates the errors, giving the results on both hard copy and punched 



48 



UJ 

oo 



O" < 

<3. h- 

o < 

CNJ Q 

•» 

X CO 

o 

n x 

V. h- 



i) 



-r 


z: 




<. 


»— i 




o 






eg 


00 






a 




X 


a: 




o 


<E 




m 


o 




s 






x 


t/> 




*» 


o 




■» 


o 




• * 


UJ 




■» 


X 




•> 


a 


m 


-t 


at: 


* 


> — « 


r* 


Q 


•> 


UJ 


UJ 



00 

i — » 

< 

H- 

C 

a 

* 









QC 


QC 


UJ 


— x 






UJ 


o 


-J 


o 






ca 




CL 


eg 






:e 


< 


s 


•***» 






ZD 


QC 


o 


H 






X 


H- 


CD 


oo 


— H 






X 




LU 


o o 




1 — 


LU 


UJ 


l— 


eg o 




UJ 




QC 




•> CNJ 




00 


UJ 


< 


•» 


CNJ 






Cf. 






II CD 




< 


< 


OO 


O 


HH K 




h— 




l — 


CNJ 


O 




< 


UJ 


UJ 




•« o 




• O 


QC 


oo 


aj 


•— < l ^ 


CD 


» 


UJ 






M — — X 


-X 


m *> 


O X 


c 


«x> 


w 0-4 


O 


X 


O Jr- 


h- 




| . 0 — » 


o 


•'in 


rO - — 


< 


o 


oOst-4-H 




<in 


•* 


O 


eg 


1— ' LU< II oo 


11 




O X-x < 


O 


— 


00 hONQUJ 


Q 


*>s 


H- nJ*» O 


•» 


< 


LU — 


2 


IZX 


cd m •• <c\j 


X 




h- C" — •'UJ • 


UJ 


— X 


o— * 


— X 




^<Q-i •'LU 


«x> 


oo 


— OOCDOX 


oo 


zz 


oor^ioz 


o 


mx 


_jino 


SOX 


u 


— 1 CNJ •'i— 1 >— 1 • 


r— 4 


•'—4 0 


• •'r-i — i *>rn 


O •'— < 


> — ( 


’ VO '--' 






0-0 — 


00 


LTU — •w-j — LO~ « 


in 


^-1 


UJ LU *x-» | — 


^'-'h 




wCLUC^w 


— 4 LU <L 


• UJC IJJC 


UJ< 


UJ 


a2h>;Q< 


O -M- *>: O 


SH 5Z0Oh- 


OH-2Z 



MU.<^*o<iZwn:h- w ^Li »H->~«r£Cja 

>-h || UJLJOCLJIUIL II UJ 1 1 r X(JUU.0Ca QT. U II OQC Oh- 
O ^ CC X IS X QZ >-« 2L (X ^ -£ X O 3: U- I— JC t JL £1 O _£ LL oo UJ 



* 


O 


OO 


CD 


o 


O 


CD CD 


oo 


-fr 


H 


ego 


m 


o 


C* 


oin 


o-u 






*»— 4 




eg 




cn 


nT 



* 

ldldo 



49 



cards. A FORTRAN Deck set-up consists of the requisite control cards, four 
parameter cards, the data deck, and finally two trail cards. 



1. TEXT Any identifying alphanumeric remarks desired (20A4) 

2. K, FRE0, FREU , DF , TAU, FEA, RH0 , THETA (110, 7F10. 5) 

K=1 Program control parameter signifying the start 

of a spectrum. 

FRE0 Upper spectrometer energy limit for the spectrum 

in MeV . 

FREU Lower spectrometer energy limit for the spectrum 

in MeV. 

DF Spectrometer energy step width in MeV. KIK will in- 

terpolate between missing data sets. 

TAU* Counting system dead time/pulse length. 

FEA* Incident accelerator energy (MeV). 

RH0* Target thickness (g/cm^). 

THETA* Electron scattering angle (degrees). 

All parameters followed by an asterisk are for identification purposes only 



and are not used in the calculations. They may be left blank if desired. 



3. a(l) . . .a(10) 



(1014) 



Spectrometer counting system normalization factors 
relative to counter number 5. The a_^ ! s are defined by: 



a . 

i 



ENERGY (5) 
ENERGY (i) 



1 * 10 4 



4. Input/Output Control 

LL, LU, L0 , LF0R, INV, FKXAD, FKEN (512,2F10.5) 

LL=1 

LU, L0 The spectrum of counter I specified by LU<I<L0 will be 
punched out onto cards. 

LF0R=O 

INV>0 The order of the input data is reversed; i.e., from 

lower to higher energies. 
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FKXAD>0 


The input capacitor voltage (full-scale setting) 
will be scaled to units of 10“^ volts by internal 
multiplication by this factor. For the NPS LINAC 
Facility, FKXAD=100 with the voltage being record- 
ed in millivolts. 


FKEN>0 


The input energy will be scaled to units of KeV 
by internal multiplication by this factor. For 
the NPS LINAC Facility, FKEN=i00 with energies 
recorded in MeV. 



The data cards now follow with two cards per data set. The actual format 
used is quite complex, since the program reads data around the colons and 
reconstructs decimal numbers. The general order of the input data fields 
is as follows: 

Card 1: E0, VFS, E, C, T, V, DUMMY, C0 (1) . . . C0 (3) 



Card 2: C0( 4) 


...C0(1O), B, A 


E0 


Incident accelerator energy 


VFS 


Full-scale voltage setting on t-he capacitor (mv). 


E 


Spectrometer energy setting (MeV). 


C 


Canacitance of the collecting canacitor QjiF). 


T 


Elapsed time (sec) . 


V 


Final voltage on the capacitor (mv) . 



DUMMY = 000000 

C0(I) • Spectrometer counter readings. 



B 


Spectrometer backing counter reading. 


A 


Spectrometer ’accidentals 1 counter reading. 


5. Trail Cards. 





When a card is encountered with columns 1-7 left blank, that card and 
the immediately following card are interpreted as control cards, with the 
following meanings: 



VFSC0 


RMIN=VFS*10-3 


=0 


RMIN=10~ 4 
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T wo points differing from each other in energy by 
less than RMIN will be included as one point on 
the spectrum. 

C0(I)>O Counter I will not be considered in compiling the 

total spectrum. 
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PLACE 2 TRAIL CARCS HERE 



C . SPECTRUM 



SPECTRUM uses the output cards from KIK as its input data, displaying 
the spectrum as a ’point plot 1 drawn on the NPS Computer Facility’s CALCOMP 
Plotter. The graph will be 15 inches long and 9 inches high. SPECTRUM re- 
quires three parameter cards; the first two cards contain alphanumeric iden- 
tifying text in columns 1-48 and the third is as follows: 

NPTS, XSCL , ESCL, E0FFST , X0FFST (13 , 7X , 4F10 . 0) 



NPTS 


Numer of points to be plotted, with the restriction that 
NPTS>30n + 1; n is any positive integer. 


XSCL 


Ordinate scale in units/inch (only one significant digit.) 


ESCL 


Abscissa scale in similar fashion. 


E0FFST 


Position of the abscissa zero point. Usually for the 
elastic resonance E0FFST=E e ^ + 2; for the inelastic spectrum 
E0FFST=E e j which results in the points being displayed at 
excitation energies. 


X0FFST 


Position of the ordinate zero point. 
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THIS PROGRAM PLCTS.OUT THE INELASTIC SPECTRUM WITH THE ENERGY SCALE CFFSET BY 
THE APPROPRIATE AMOUNT. THE ELASTIC PEAK WILL USUALLY EE LCCATEC AT THE 
ORIGIN SO THAT THE INELASTIC PEAKS MAY BE READ DIRECTLY IN EXCITATION ENERGY. 
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SAMPLE FGRTR AN CECK SET UP 
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D. NAW 



NAW uses the output from KIK as its input data. This program calculates 
the ratio of the area under the inelastic resonance peaks to the area under 
the elastic resonance peak. Because of its extreme length (3134 FORTRAN 
cards) the program listing will not be reproduced; documentation may be 
found in reference 15. NAW accepts up to 198 elastic data points, 498 in- 
elastic data points, and 25 inelastic resonances per spectrum. 



1. Text Alphanumeric identifying remarks. (20A4) 

2. EICH, PN 2 , PN 3 , PN 4 , PN 5 , PN 6 , PN ? (7F10.5) 

EICH Spectrometer energy normalization factor to scale all 
input energies to MeV. 

PN. Starting values for the various parameters being used to 

fit the elastic spectrum with an analytic expression. They 
are used only when L^>0 on card //7. 

PN 3 Half-width of a gaussian to fit the high energy side of 
the elastic peak. 

PN^ Half-width of a gaussian to fit the low energy side of the 
elastic peak. 

PN 4 Matching point for the low energy gaussian and a hyperbolic 
segment. PN^ is specified in MeV from the position of the 
zero point and furthermore is a negative number. 

.. PN^ Position of the zero point. When PN^=0 the zero point is 

taken to be the position of the elastic resonance peak. 

PNg>0 Height of the second elastic resonance peak. 

PNy>0 Position of the second elastic resonance peak. 

3. EE, THETA, D, UNT, GINT (5F10.50) 

EE Elastic resonance energy in MeV. 

THETA Electron scattering angle in degrees. 

2 

D Target thickness in g/cm~. 

UNT Constant background at the high energy end of the elastic 
resonance spectrum in counts //iC. 
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GINT 



Lower integration limit for calculating the area under the 
elastic peak. This value is only used when L^=0 on card // 7. 

At all other times the lower limit is fixed at E e ^ minus twice 
the peak width. 

4. EEL , EEL , GSSW , GSSW , DSSW, DWAL (6F10.5) 

o y U y o y u ’ y v ' 

EEL Q High energy limit for fitting the elastic peak. 

EEL U Low energy limit for fitting the elastic peak. 

GSSW Q Upper limit for calculating the radiation tail. 

GSSW U Lower limit for calculating the radiation tail. 

DSSW Energy step to be used in calculating the radiation tail. 

DWAL See 1>2 on card //7. 

5. E q , E u (2F10.5) 

E q High energy limit for fitting the inelastic spectrum. 

E u Low energy limit for fitting the inelastic spectrum. 

6. GRl, GRl, GRg, GR2, Gr3, GR3 (6F10.5) 

Regions where it is expected that only background contributes to 
the spectrum; i.e., between resonances in the inelastic spectrum. 
There must be at least one set GR Q , GR U input, with a maximum of 
three such sets. 

7. L-^ , li2 > . L 22 (2213) 

L^ Analytic expression used to fit the elastic resonance. • 

=0 Simple curve sections (ARC3) . 

=1 A continuous curve consisting of two gaussians and a seg- 
ment of a hyperbola (ARC5) . 

=2 The elastic resonance will be fit as if Lj=l and the in- 
elastic resonances will be fit using the Breit-Wigner formula. 
Each curve will be multiplied by a straight line with equa- 
tion Y=1 + Pi(Ex~E ox )* Hi will be specified on card //10; 

E x is the excitation energy, and E QX is the position of the 
maximum. 

=3 The elastic resonance will be fit as if L-^=2 except that 
NAW will determine the values and card //10 will be left 
blank. 

L 2 Analytic expression used to represent the background due 
to the radiation tail. 
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=1 BG = Pj^SSW 

=2 BG = P].*SSW + P 2 

=3 BG = P X *SSW + P 2 + P 3 X 

=4 BG = P x + SSW 

=5 BG = P x + P 2 X + SSW 

=6 BG = P x + P 2 X + P 3 X 2 + SSW 

=7 BG = P x + P 2 A + SSW 

=8 60=?!+ P 2 A + P 3 A 2 + SSW 

=9 BG = P x + P 2 B + P 3 B 2 + P 4 X + SSW 

=10 BG = Px + SSW*(P 2 + P 3 X) 

= 11 BG = Pj_ + SSW*(P 2 + P 3 X + P 4 X2) 

will be determined by NAW. 

SSW = Theoretically calculated radiation tail. 

X = E - Eg. The reference energy Eg will be internally calcu- 
lated by NAW. 

A = l-(E el “ E). 

B = 1/ (E e i - E - DWAL) . Should L 2=9 be selected, then DWAL must 
be input on card //4. If Lg^9 then the space for DWAL on card // 4 
should be left blank. 

Lg The number of resonances with specified energies. If L^<0 
then the resonances will be kept fixed relative to E x (l) , 
but E x (l) will be positioned by NAW. 

The number of resonances with energies to be determined by 
NAW. Note that Lg + must not exceed 25. 

Lg The number of line widths to be determined by NAW. 

=0 All line widths will be specified on card //9. 

>0 E.g. =3, means that three line widths are free and in fact 
are E x (l), E x (2), and E x (3). If there are more than three 
resonances, *E X (3) through E X (N) will all be fitted with the 
same line width. 

L 6 =0 
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L -j Target geometry: 

=0 Transmission - normal thickness entered on card // 3. 

=1 Reflection - normal thickness entered on card //3. 

1 Effective thickness entered on card //3. 

Lg Output control: 

=0 Maximum output. 

=1 Minimal output - no tables and no intermediate calculations. 
=-l Limited output - no tables. 

Lg Plotting control: 

=0 Results are not plotted. 

=2 Results are plotted. 

Lfo Repetitions of the NAW procedure. 



=0 No repetitions are performed. 



>0 In this case it is possible to evaluate the data with new 

parameters. The new parameters will be entered on additional 
cards placed after the data. 

^ll~Li4 Inelastic spectrum fits iteration control. Only override 
to the default options need be entered. 



L ll=2 

L 12 

L 13 

L 14= 6 

Li5=° 

<0 

L 16“ L 18 

L 16 



=1 

=2 



Maximum number of iterations; default is L^2=30. 
Minimum number of iterations; default is L^-5. 



The radiation tail will be calculated using Ginsberg* s method. 

The radiation tail will be calculated using Schiff*s method. 

Radiation tail fit control. Only override options need be 
entered. 



Analytic representation of the theoretically calculated 
radiation tail; default is L^=l. 

SSW = P, + P 0 A + P q A 2 + P,A 3 + P,A 4 
1 2 3 4 5 

SSW = + P 2 A + P 3 A 2 + P 4 E + P 5 E 2 
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=3 SSW = ? 1 + P 2 A + P 3 E + P 4 E2 + P 5 e 3 
=4 SSW = PjA + P 2 E + P 3 E 2 + P 4 A 2 + P 5 E 3 
=5 SSW = P 2 + P 2 U + P 3 U 2 + P 4 U3 + P 5 U^ 

P^ will be determined internally by NAW. 

A = l/(E el - E) 

U » E - Eh 

E b = GSSW U + (GSSW q - GSSW u )/L lg 

<0 The radiation tail parameters and E^ will be specified 
on card #15. 

L 17 The number of parameters for the radiation tail fit 
( 3 ~Li 7 - 5 ); default is Lj 7 = 4 . 

Lfs See default is L-^g=2. 

L^g>0 Error ellipsoid control values will not be calculated near 
the fit. 

l 20 =0 

L, 2 i >0 The calculation of the radiation tail will be modified by 
the factors entered on card #16. 

L^^>0 Inelastic spectrum results will be punched out onto cards. 

8 . E x (1) > E x (2) > •••» E x (7) ’ LE (7F10. 5,110) 

Resonance excitation energies in MeV ; enter those with 
variable .position first. For more than seven resonances 
set LEX) and then continue listing the additional reso- 
nances on successive cards. 

8 a. E x ( 8 ) ... E x (14), LE 

8 b. E x (15) ... E x (21), LE 

8 c. E (22) ... E (25) 

A A 

9. V(l), V(2), ..., V(7), LBR (7F10.5,I10) 

Resonance line widths. For more than seven resonances set 
LBR>0 and continue on successive cards as on card # 8 . If 
some or all V(I) are zero, then the line width used will be 
the dispersion of the spectrometer counting system. 

10. P(l), P(2) , ..., P (7 ) , LP (7F10 . 5 , HO) 

The P^ factors used when L^=2 on card #7. Set LP>0 for 
more than seven factors and continue as for card # 8 . 
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11. Z, A, XRAD, RMS, R, G, C, T 



(8F10.5) 



Z Atomic number. 

A Mass number. 

XRAD Radiation length in g/cm^ 

RMS Root mean square nuclear radius in fm. Used only when N=2 
on card #14, otherwise left blank. 

R, G When N=0 or N=1 on card #14, R and G are interpreted as 

the parameters in the Helm expression for the form factor. 
When N=3 , R=CRT and G=TRT . When N=2 or N=4, R and G are 
left blank. 

C, T The fermi charge distribution parameters. Since these val- 
ues are used as the parameters for the trapezoidal charge 
distribution potential used in the radiation tail calcula- 
tions, the values must be entered even when N/4. 



12. XMAS, H 

XMAS Nuclear mass equivalent in MeV. 

H Isotope percentage (100% = 1.00) 

13. A, B, C, D, E 



(2F10.5) 



(5F10.5) 



The A ... E are the coefficients for the calculation of; the 
elastic form factor using a power series expansion onq : 



■el 
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When A=0 the elastic form factor will be calculated using 
the Born Approximation. 



14. A, N 

A • Mass number. 



(I3,24X,13) 



N Selection of the nuclear model to be used in the Born 

Approximation to calculate the form factor: 

=0,1 F = 3(j 1 (qR)/qR)exp(-G 2 q 2 /2) Helm. 

-q 2 <R 2 >/6 + q 4 <R 2 > 2 /120 + q 6 <R 2 > 3 /5040 Reihe. 

_ nD'T*.2\„„„ / 'T'DT*„2\ nl.ll 



(6F10.5) 



=0, 


1 F = 


=2 


F = 


=3 


F = 


=4 


The 


PI, 


P2, 



This card is read only if Ljg<0. 
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P1-P5 Radiation tail parameters, see Lj^ for the expressions. 

This value is read only when Ljg=-5. 

16. XSSW, XM, XB, LSSW (3F10.5,I3) 

This card is read only if L^X). 

XSSW Factor multiplying the area under the elastic peak. 

XM Factor multiplying the Holler term. 

XB Factor multiplying the Bremsstrahlung term. 

The equation for the radiation tail is now given by: 

(Max I) 

SSW = ( — ^ “ + XM(Moller term) + XB(Brehm. term) ) • (XSSW)A £ 

F a M 

LSSW>0 The calculated radiation tail will be punched onto cards. 

The data cards now follow. The elastic and inelastic data must 
be separated by a card having the data energy step width entered 
beginning in column 72. Repeated evaluations of the same set of 
data (i.e., L^q> 0) must have the following set of cards placed 
after the data cards - one set for each repetition to be performed 



1) 


Lj ... 


L 10 


see 


card 


n 


2) 


E c E „ 




see 


card 


#5 


3) 


MD • 


•• E X (N) 


see 


card 


#8 


4) 


v(l) .. 


. V (N) 


see 


card 


it 9 


5) 


P(l) .. 


. P(N) 


see 


card 


no 
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E. MARPHASE 



MARPHASE uses the Rawitscher-Fischer phase shift code YALERF developed 
at YALE [Ref. 8j to calculate the elastic electron scattering cross sections, 
Mott cross sections, and form factors. A parameter set consists of two 
cards. 

1. Z, A, PI, P2, P3, MODEL (5F10. 5,110) 

Z Atomic number. 

A Mass number. 

P3=0 

M0DEL=1 Fermi model of the nucleus is assumed with 

P(r)= : — !— — ■ — . 

I ♦ exp [ r - c/ct/4.4.)) 

In this case Pl=c and P2=t, both in fermis. 

=2 The Shell model of the nucleus is assumed with 

p(r)= [ + a(xf] ^P[-(a)) ‘ 

In this case Pl= a and P2=A. 

2. E0, Tl, T2 , TH, R0, LAG (5F10.5) 

E0 Electron beam incident energy in MeV. 

Tl Starting scattering angle in degrees. 

T2 Final scattering angle in degrees. 

TH Angular increment in degrees. 

R0 The fitting radius pertaining to matching the Nuclear 
and the Coulomb wave functions. 

LAG The number of phase shifts. 

3. Another card #2 placed next will cause the program to use the 
same information for card #1 to compute the new combination of incident 
energy and scattering angles. A blank card is interpreted to mean that 
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a computation is now desired for a different set of card # 1 parameters and 
is followed by a new parameter set, usw. 
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F. DUKE 



DUKE calculates the theoretical cross sections for both electric and 
magnetic giant resonance transitions. Either the DWBA or the PWBA may be 
selected. A parameter set consists of four cards. A subsequent parameter 
set consists of only the last three cards. A blank card signifies the end 
of the desired calculations. 

1. Z, AX, Zl, El, E2 , ETRAN, B0RN, M0RE (7F10.5,I10) 

Z Atomic number. 

AX Parameter c for the Fermi charge distribution. 

Zl Parameter t for the Fermi charge distribution. 

El Starting energy in MeV. 

E2 Final energy in MeV. 

ETRAN Excitation energy in MeV. Set ETRAN=0 for a pure coulomb 
transition and set ETRAN=E1-E2 if the transverse component 
of the cross section should be considered. 

B0RN=O Distorted Wave Born Approximation. 

=1 Plane Wave Born Approximation. 

M0RE=O Electric transition. 

=1 Magnetic transition. 

2. L, KMAX, RADVAR, SKIVAR, ASG , ZKMAX, ETRINT, LM0D (2110 ,5F10 . 5 , 110) 

L Multipolarity of the transition. 

KMAX>0 Number of partial waves. 

=0 KMAX is calculated internally: KMAX =9+ (El-30) /10. 

RADVAR Position of the transition distribution normalized to c. 
SKIVAR Width of the transition distribution normalized to t. 

ASG=0 Complete output is printed. 

=1,2,3 Limited output is printed. 
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ZKMAX This modifies KMAX, as computed internally, to be 
KM AX = KMAX + ZKMAX. 

ETRINT Parameter used to determine the integration steps. ETRINT = 

35.0 by default; increasing ETRINT decreases the step width, 
increases the accuracy and increases the execution time. 

LM0D=O Transition distribution function will be approximated by 

the derivative of the ground state fermi distribution func- 
tion. 

=1 Transition distribution function will be approximated by the 
ground state fermi distribution function. 

3. HHl , HII2 , NST1 , NST2 , NST3 , IST1 (2F10 . 5 ,4110) 

This card is best given as a blank card, DUKE then determines the 
integration parameters using the internal SUBPROGRAM INDEX; see 
below for assumed values. 

HHl Integration stepwidth inside the potential in fermis. 

HH2 Integration stepwidth outside the potential in dimensionless 
units . 

NST1 Length of the internal region in units of AX. 

NST2 Length of the external region in dimensionless units over 

ETRINT. 

NST3 The values of the radial wave functions will be computed in 

three regions delimited by: NST1, HHl + NST2 , HH2 + i’NST3, HH2; 

where i = 1, 2, ..., IST1. 

When this card is left blank then: 

HHl = .25 fm. 

HH2 = 0. 3t)c/ ^ El • E2 

NST1 = 3AX/HH1 

NST2 = (30. + (ETRAN/E1)- ETRINT) /HH2 

NST3) Mean of the integrated radial wave function within one wave- 
IST1J length of the oscillation with the largest wavelength. IST1 = 6, 
and NST3 = wavelength/ (HH2 • IST1) 

4. THEMN, DELTH, N0, STANZ, NZDG (2F10. 5 , 1110 ,F10 . 5110) 

The cross section for a parameter set will be calculated for all 
angles 0 = THEMN + i-DELTH, i=2 , 1, . . . ,N0. KMF, STANZ, and KZDG 
are test parameters and have no effect when left blank. If this 
card is left blank then THEMN = 30, DELTH - 15, and N0 = 10. 
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G. BCALC 



BCALC is used to calculate the transition strength, B. While written 

specifically for use with 75 MeV and 90 MeV incident energies, it is easily 

changed to accommodate other values. BCALC uses the following data cards: 

one card with the number of 75 MeV experimental points; a card for each such 

point containing the scattering angle in degrees and the experimental, value 
2 

of F , FORMAT (13 , 7X,E10. 4) ; a similar set-up for the 90 MeV experimental 

2 

points; and lastly, the theoretical values for scattering angle and F . 

For each transition type to be evaluated a set of theoretical values must 
be supplied - 75 MeV values followed by 90 MeV values. The output from 
BCALC is a listing of all experimental values with the corresponding the- 
oretical values printed next to them, the computed value for B, the com- 
puted value of the quantity VAR (see Section IV-A) , and finally a listing 
of the theoretical value of F multiplied by B. 
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work, measurements were made to identify the giant nultipolc resonances of 
^9 ^Au. Five inelastic scattering spectra wore obtained and analyzed for 
beam energies and angles of 75 McV / 75°, 75 MeV / 90°, 75 McV / 120°, 

90 MeV / 75°, and 90 McV / 90°. In each spectrum, giant muJtipolc resonances 
wn i ■ I a l ■ yi ir.ii I r> ■ »■ *r >*i . .1 '+1, h . 0 and 14 McV. The 

strengths and multi polarities were found to he B(M3 ) - 1 3 , OOOu 7 fin 1 , 
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B(E2)=8,000e 2 fm 4 , B (E2)=9 , OOOe 2 fm 4 , and B (El)=90e 2 fm 2 , respectively. The 
full widths at half maximum for these resonances were found to be, respect- 
ively, 3.1, 2.9, 3.1, and 4.9 MeV. 

The resonances found at 10.9 and 14 MeV confirm the only other measure- 
ments made on gold £ 4 ]. The resonance found at 7.1 MeV agrees with those 
findings concerning the magnetic character of the transition but differs 
in the multipolarity assignment; however, our assignment may also be con- 
sidered tentative since insufficient data were obtained at the backward 
scattering angles. The resonance found at 9.1 MeV has not been previously 
reported . 
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